Locking electron spins into magnetic resonance by electron-nuclear feedback 



I. T. Vink, K. C. Nowack, F. H. L. Koppens, J. Danon, Yu. V. Nazarov, and L. M. K. Vandersypen 
Kavli Institute of Nanoscience, Delft University of Technology, 
PO Box 5046, 2600 GA Delft, The Netherlands 
(Dated: February 16, 2009) 



The main obstacle to coherent control of two- 
level quantum systems is their coupling to an un- 
controlled environment [1]. For electron spins 
in III-V quantum dots, the random environ- 
ment is mostly given by the nuclear spins in 
the quantum dot host material; they collectively 
act on the electron spin through the hyperfine 
interaction, much like a random magnetic field 
[2, 3, 4, 5, 6, 7, 8]. Here we show that the same 
hyperfine interaction can be harnessed such that 
partial control of the normally uncontrolled envi- 
ronment becomes possible. In particular, we ob- 
serve that the electron spin resonance frequency 
remains locked to the frequency of an applied 
microwave magnetic field, even when the exter- 
nal magnetic field or the excitation frequency are 
changed. The nuclear field thereby adjusts itself 
such that the electron spin resonance condition 
remains satisfied. General theoretical arguments 
indicate that this spin resonance locking is accom- 
panied by a significant reduction of the random- 
ness in the nuclear field. 

In thermodynamic equilibrium, the nuclear spins in 
the quantum dot host material are randomly oriented, 
even at dilution refrigerator temperatures and in mag- 
netic fields of a few Tesla. An electron spin confined 
in the quantum dot interacts via the hyperfine coupling 
with A'' ~ 10^ nuclear spins and as a result experiences 
a random nuclear field Bjq. This random nuclear field 
is sampled from a distribution with a root mean square 
width oc A/ g^B"/N ^ where g is the electron ^-factor, /is 
the Bohr magneton and A the hyperfine coupling con- 
stant [k, 135/xeV in GaAs). Measurements typically give 
a width of ^ 1 mT. As a result, we lose track of the phase 
of a freely evolving electron spin within a time Tj* of a few 
tens of nanoseconds i, II IS 0, i- Similarly, when the 
spin evolves under an oscillating driving field, the nuclear 
field leads to a random offset in the resonance condition 
which has a comparable amplitude to presently achiev- 
able driving fields. This results in poorly controlled spin 
rotations [9|]. 

It is therefore of great importance to develop the abil- 
ity to control and manipulate the nuclear field with great 
precision. In particular, it would be highly desirable to 
set the nuclear field to a narrow distribution of values at 
the start of every experiment 0, El, H m • This would 
immediately reduce the rapid dephasing, and the electron 
spin would loose phase coherence only from the slow sub- 



sequent evolution of the nuclear fiel d, g iving a predicted 
spin coherence time of 1 — 10/j,s lj,ll5|. Such narrowing 
has been achieved in an ensemble of self-assembled quan- 
tum dots by synchronizing the precessing spins with a se- 
ries of laser pulses [3| • More recently, the spread of the 
difference in nuclear fields in two neighbouring quantum 
dots was reduced via a gate voltage controlled pumping 
cycle, giving a 70-fold increase in the Tj* for states in the 
two-electron = subspace [l3| . 

Here we exploit electron-nuclear feedback in order to 
control and manipulate the nuclear fields in two coupled 
quantum dots during continuous wave (CW) driving of 
the electron spins in the dots. We observe that the nu- 
clear field adjusts itself such that the electron spins re- 
main in resonance with a fixed driving frequency even 
when we sweep the external magnetic field away from 
the nominal resonance condition. Similarly, the electron 
spin resonance frequency remains locked to the excita- 
tion frequency when the excitation frequency is swept 
back and forth. These distinctive features set our ob- 
servations apart from the many previous observations 
of dynamic nuclear spin polarization in quantum dots, 
both in transport ^ , Il9l . |20| . 21 , and optical mea- 
surements [U 0, llal ■ We investigate the origin of this 
feedback by studying its dependence on the amplitudes 
of the applied ac magnetic and electric fields and on the 
sweep rates. We show theoretically that the spin reso- 
nance locking must be accompanied by a narrowing of 
the nuclear field distribution, in the present experiment 
by more than a factor of 10. 

The measurements are performed on an electrostati- 
cally defined double quantum dot tuned to the Pauli spin 
blockade regime [l^l , with effectively one excess electron 
on each dot (the actual electron number is small but un- 
known). When the two electrons have parallel spins, the 
electron flow through the dots is blocked. When one of 
the spins is flipped, the spin blockade is lifted and elec- 
trons flow through the two dots until the system returns 
to a state with parallel spins on the two dots. As previ- 
ously demonstrated Q , it is possible to flip the dot spins 
via magnetic resonance, by ac excitation of an on-chip 
wire which generates an oscillating magnetic field at the 
dots: when the excitation frequency, /, matches the elec- 
tron spin resonance (ESR) frequency, \g\^BBf)/h, a finite 
current fiows through the device. Here h is Planck's con- 
stant, and _Bo the external magnetic field. In addition, 
current can flow at zero magnetic field, where the elec- 
tron spins can flip-flop with the nuclear spins in the sub- 
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strate [l9[. We use this zero-field feature to determine 
and adjust for small magnetic field offsets present in our 
setup. The zero- field peak and the ESR response are 
seen in current measurements under CW excitation with 
increasing excitation frequency at fixed magnetic fields 
(Fig. [T^, similar to the data pubhshed in Ref. Q , and 
taken on the same device but in a different cooldown. 

Surprisingly, when we reverse the sweep direction, a 
distinctly different behavior is observed over a wide range 
of dot settings (see Supplementary Information for details 
of the tuning parameters). Current starts flowing when 
the driving frequency hits the spin resonance frequency 
but remains high even as the frequency is swept well be- 
low the nominal resonance condition (Fig. [Tja. The fact 
that the current remains high implies that the electron 
spin is still on resonance with the excitation frequency, 
and that an effective field, Seff, counteracts the exter- 
nal magnetic field Bq: hf = \g\nB{Bo + Bes). From the 
fact that the current is strongly reduced when we simul- 
taneously excite any of the three nuclear spin species in 
the substrate (data not shown), we conclude that this 
effective field is created by dynamical nuclear spin polar- 
ization. This nuclear field builds up exactly at the right 
rate in order to keep the electron spin in resonance with 
the changing driving frequency, which implies there is a 
built-in electron-nuclear feedback mechanism. 

Similar dragging of the resonance is observed when 
sweeping the magnetic field for a fixed excitation fre- 
quency. In Fig. [2^ the magnetic field is swept from -33 
mT to 97 mT (right vertical axis) in about 25 seconds. 
We first see the zero-field peak, as expected, and next the 
current jumps up around Bq — 67 mT, which is slightly 
below the nominal resonance condition (/ ~ 400 MHz, 
\g\ — 0.36). The current remains high as the field is 
swept further to 97 mT, which is well outside the ESR 
linewidth in the absence of feedback (see Fig. |4)d). Sim- 




ldo,(PA) 




-100 -50 50 
Magnetic field (mT) 



-50 50 
IVlagnetic field (mT) 



FIG. 1: Electron spin resonance locking during frequency 
sweeps, a, Current through the double dot (colorscale) sub- 
ject to CW magnetic excitation, when sweeping the frequency 
up at fixed magnetic fields. The bright fork indicates the po- 
sition of the ESR condition, b, Similar to a but sweeping the 
frequency down. The ESR frequency remains locked to the 
excitation frequency when the excitation frequency is swept 
past the nominal resonance condition. 
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FIG. 2: Electron spin resonance locking during magnetic field 
sweeps, a. Current through the double dot as a function of 
time, while the magnetic field is first ramped up (right axis) 
and subsequently held fixed, under CW excitation (/ = 400 
MHz), b-c. Two current traces similar to a, but after the 
magnetic field is ramped up, it is repeatedly swept down and 
back up over a 30 mT range (right axis). After the ESR 
condition is first met, the electron spin remains locked into 
magnetic resonance for up to two minutes, even though the 
resonance condition is shifted back and forth. 



ilar to the case of the frequency sweeps, a nuclear field 
builds up exactly in such a way as to maintain the ESR 
frequency locked to the excitation frequency. When we 
subsequently keep the field fixed at 97 mT, we observe 
that the electron spin can remain locked into magnetic 
resonance for well over a minute. 

It is also possible to drag the nuclear field back and 
forth under fixed-frequency excitation. In Figs, ^jp and 
Bq is ramped up from -33 mT to 117 mT, and is sub- 
sequently swept back and forth between 117 mT and 87 
mT in a triangular pattern. The current again jumps up 
as we sweep through resonance and subsequently remains 
high independent of the sweep direction. In Fig. [2}: the 
resonance is lost after approximately 1 minute, whereas 
in Fig. the spin remains locked on resonance during 
the entire experiment (about 2 minutes). 

These remarkable observations of spin resonance lock- 
ing due to electron-nuclear feedback are characterized by 
a number of common features. First, the current jumps 
up abruptly, in many cases in less than a few 100 ms, at 
a field value that varies over 10-30 mT around the nom- 
inal resonance condition (see the green circles in Fig. [4] 
below). This is a further indication that the system is 
actively pulled into resonance - without feedback a cur- 
rent peak with smooth flanks and a width of a few mT is 
expected 27[. Second, the resonance dragging generally 
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occurs only for fields larger than the nominal resonant 
field, or for frequencies lower than the nominal resonance 
frequency. This is opposite to the case of the usual Over- 
hauser effect, as discussed further below. Third, the ini- 
tial current jump is usually followed by a second current 
jump, before the current drops back to zero. A possible 
explanation for this double step is that the first current 
plateau corresponds to a situation where both electrons 
are on-resonance, and that only one electron remains on 
resonance after the second jump (see Supplementary In- 
formation for a discussion of the current levels) . When 
the resonance is lost in this last dot too, the current re- 
turns to zero. 

This interpretation of the double current step is sup- 
ported by pump-probe measurements shown in Fig. [31 
Starting from the second current plateau with Bq = 80 
mT and / = 276 MHz, we switch off the CW excitation 
and probe the position of the ESR frequency as the nu- 
clear field returns to equilibrium (we use short bursts 
for probing in order to minimize feedback during the 
probe phase). We see that the ESR frequency returns 
to its nominal value, slightly above 400 MHz, within 20 
seconds, corresponding to the relaxation time of the lo- 
cal nuclear spin polarization. This signal must originate 
from the dot that is still locked into magnetic resonance 
at the end of the pump phase. In addition, we see a re- 
sponse at the nominal resonance frequency already from 
the start of the probe phase. Presumably, this signal 
arises from the other dot, where the resonance was lost 
and the nuclear field has (nearly) relaxed by the time the 
probe phase starts. 

In order to better understand the locking mechanism, 
we study how far the nuclear spin polarization can be 
dragged by performing magnetic field sweeps as a func- 
tion of the applied microwave power, the microwave fre- 
quency and the magnetic field sweep rate. Specifically, 
we repeatedly ramp the magnetic field from -28 mT up- 
wards and record (i) the field at which the current jumps 
up (circle in Fig. |4^), (ii) the field where the current 
jumps to a still higher value (diamond symbol), and (iii) 
the field where the current drops back to zero (cross). 
The resulting data points are shown as scatter plots in 
Figs. Hh-e, using the same symbols. 

The first current jump always occurs as the nominal 
resonant field (in the absence of feedback) is first ap- 
proached. The second jump and the current drop occur 
at fields that increase with driving amplitude over the 
range that we could explore (for still stronger driving, 
spin blockade was lifted by photon assisted tunneling, 
and we lost sensitivity to spin flips). For the highest pow- 
ers accessible in the experiment, the electron spin is main- 
tained on resonance over a magnetic field range of a few 
100 mT. As the power is reduced, the locking effect van- 
ishes. Furthermore, the field that can be reached before 
the resonance is lost, increases with excitation frequency. 
Earlier measurements on the same sample showed that 
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FIG. 3: Pump-probe measurement of the relaxation of the nu- 
clear spin polarization. At a fixed magnetic field of Bo ~ 80 
mT, we apply CW excitation {P = —13 dBm) sweeping the 
frequency from 500 MHz to 276 MHz at 43 MHz/s, and drag- 
ging the nuclear field along (pump phase). Next we turn off 
the CW excitation and apply 140 ns microwave bursts every 
2 fis at frequency /probo throughout a 40 s probe phase. This 
pump-probe cycle is repeated for different probe frequencies, 
277 MHz < /probo < 450 MHz (see vertical axis). The hor- 
izontal axis indicates the time t into the probe phase; the 
data for t < correspond to the pump phase. In the pump 
phase, the current (plotted in colorscale) jumps up twice, 
reaching the highest current plateau (traces where the res- 
onance is lost by the end of the pump phase are left out). 
When the frequency is switched to /probe a,t t — 0, the cur- 
rent drops to zero since the excitation is now off-resonance. As 
the nuclear spin polarization relaxes, the resonance condition 
|(7|/ifl(i3o -I- -Bjv(t)) = ^/probc will be fulfilled at some point 
in time at which the current sets on again. Varying /probe 
reveals then the nuclear spin relaxation as indicated by the 
white dashed line (guide to the eye) marking the onset of the 
current. Even though the excitation is applied only in bursts, 
the electron spin nevertheless remains locked into resonance 
in some cases, stalling the nuclear spin relaxation. The orange 
dashed line marks an additional signal at the nominal reso- 
nance frequency already present from the start of the probe 
phase. 



along with the ac magnetic field an ac electric field is 
generated whose amplitude for a fixed power (and mag- 
netic field amplitude) increases roughly linearly with the 
excitation frequency The dependence on driving fre- 
quency can therefore also be interpreted as stronger lock- 
ing for higher electric field amplitudes. Finally, we see 
that for higher magnetic field sweep rates the resonance 
is lost at lower fields. 

A few basic considerations give insight in the mech- 
anism behind these observations. For clarity, we here 
present a single dot picture; the results for two cou- 
pled dots are qualitatively similar [2§|. We define x as 
the dimensionless nuclear spin polarization in the dot, 
with —1 < a; < 1 (in our experiments, <C 1). In 
the absence of any excitation, the polarization naturally 
relaxes to zero at a rate a;/r„, due to nuclear spin dif- 
fusion. However, the nuclear spin dynamics will be af- 
fected through hyperfine-mediated electron-nuclear flip- 
flops when the electron spins are brought out of equi- 
librium 18, [l^ [23|. In the spin blockade regime at 
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FIG. 4; ESR locking dependence on excitation power, fre- 
quency and sweep rate, a, Current through the double dot as 
the magnetic field is swept up (/ = 400 MHz), b, Similar to a 
but now Bo is swept down. No dragging eff'ects are observed; 
the narrow peak gives the position of the nominal resonant 
field, c, Scatter plot of the switching fields as indicated by 
the symbols in a, as a function of the power applied to the 
on-chip wire. The corresponding resonant magnetic field am- 
plitude Bi at the dot is given as well, d, Scatter plot similar 
to c, as a function of /. The estimated electric field amplitude 
El increases with /, and is shown in the figure, e, Scatter plot 
similar to c as a function of magnetic field sweep rate. Blue 
lines: average and standard deviation of the magnetic fields 
where the second current jump is observed. Purple curve: fit 
of these average values with a theoretical model (see Supple- 
mentary Information). We note that there is no build-up of 
Bn in the limit of zero sweep rate, so the predicted switching 
field first increases with sweep rate, before decreasing. 



finite Bq, such non-equilibrium dynamics occurs when 
the electron spins are resonantly excited by an external 
microwave magnetic or electric field. Regardless of the 
relevant microscopic processes, we thus expect in very 
general terms a polarization-dependent pump rate Tp{x), 
which is non-zero only close to the resonance condition 
Ax'''^ = gtisB'^'' = \9\i-).bBq - hf. Here is thus the 
nuclear spin polarization that brings the electron spin 
in resonance with the excitation. The dynamics of the 
polarization in the dot is then described by 



(1) 



Fig.[5]qualitatively visualizes Eq.[l]in the form of a pump- 
ing curve for three different values of x'-"^ , where we have 
(for now arbitrarily) chosen the resonant contribution 
to be positive. From the figure we can see that stable 
points of nuclear polarization occur when dx/dt crosses 



zero with a negative slope: if x is higher (lower) than the 
stable polarization xq, dx/dt is negative (positive) and x 
gets pushed back to xq. Due to nuclear spin relaxation, 
there is almost always a stable point at a; = 0. Depend- 
ing on the particular shape of rp(a;), hence on the specific 
experimental regime, there can be one or more additional 
stable points [2i, 0, [s^l • 

We now interpret the field sweep experiments within 
this simple picture. First, given that the current remains 
high in field sweeps, a stable point must exist close to res- 
onance, in agreement with our expectation of a resonant 
peak in Tp. Next, since dragging is generally observed 
only for x > 0, Tp must be positive, as in Fig. [Sj Fi- 
nally, from the maximum nuclear field B'^'^^ that can be 
achieved by dragging, we can estimate the height of Tp-. 
when the maximum of the pumping peak falls below zero, 
i.e. when nuclear spin relaxation dominates the resonant 
pumping, the stable point at x > disappears and Bpf 
relaxes to zero (Fig. [5l red curve). 

During actual field sweeps, the resonance is lost at 
fields below B'^^^: since a dynamic equilibrium is reached 
when dx/dt = \g\^BBQ/ A instead of dx/dt = 0, the sta- 
ble operating point moves up the pumping curve (see 
Fig. O and disappears when the sweep rate exceeds the 
maximum of the pumping peak. In practice we will loose 
the resonance even earlier, because intrinsic nuclear field 
fiuctuations can drive the nuclear field across the maxi- 
mum. We model the average switching field taking into 
account such fluctuations by assuming an exponential de- 
pendence of the switching rate on the "barrier height" . 
The result is illustrated in Fig.|4^. This combined picture 
captures very well the experimental observation that for 
higher sweep rates the resonance is more easily lost, but 
not at exactly the same field every time. 

We next turn to the nature of the extrinsic pumping 
process, Tp. First, the stable points in the experiment 
generally occur for a; > 0, i.e. the nuclear field points 




FIG. 5: Nuclear spin pumping curves. The nuclear spin polar- 
ization rate for one dot [dx/ dt) is shown as a function of its po- 
larization X. The overall negative slope is due to nuclear spin 
relaxation and the resonant peak is due to the external driving 
(the green, blue and red curves correspond to three different 
^rcs-j Qij-Qies indicate stable points in nuclear spin polariza- 
tion and are found whenever the curve crosses the x-axis with 
a negative slope. During a field (or frequency) sweep, a dy- 
namic equilibrium is reached where dx/dt = \g\^BBQ/A. 
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against the external magnetic field. This is opposite to 
the usual Overhauser effect, where electron spins are ex- 
cited by magnetic resonance and relax back from j to f 
by flip-flopping with the nuclear spins. "Reverse" pump- 
ing is possible when there is an excess of t electrons, 
which are excited to | by resonant electric fields, whe reby 
the nuclear spins absorb the angular momentum 3^ 31 1 . 
In our experiment, spin-exchange with the leads due to 
photon-assisted tunneling (at 10 — 100 kHz) gives an ex- 
cess of t electrons, which favors reverse pumping. Sec- 
ond, the locking effect gets stronger, hence Tp becomes 
larger, not only with stronger driving in general (Fig. [41:), 
but also with stronger electric excitation by itself (higher 
/, Fig. mi). Based on these observations, we suggest that 
electric-field assisted electron-nuclear flip-flops combined 
with electron spin relaxation are mainly responsible for 
the resonant pumping [281] . 

Finally, we analyze theoretically the implications of 
our observations for the width of the nuclear field dis- 
tribution. We define r±(x) as the total positive and 
negative nuclear spin fiip rates that result from the in- 
trinsic relaxation and resonant response combined, so 
dx/dt = ;^(r_|_ — r_). We also define 7(2;) as the to- 
tal rate of nuclear spin flips, 7 = ;^(r+ + r_). Using 
the fact that the pumping curve exhibits a resonant peak 
at \xo\ ^ 1, we can then approximate the variance of 
the nuclear polarization distribution around xq as (see 
Supplementary Information) 
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The numerator is the local diffusion rate, and the de- 
nominator is the restoring force - the steeper the slope 
of dx/dt, the stronger the restoring force. For the case 
without pumping, we have r± = 7Vzp/2T„, so Eq. [2] gives 
us the usual result ct^ — 1/N. For a stable point xq > 
near resonance, we take as a rough estimate for the lo- 
cal slope the maximum of Tp divided by its width. This 
gives ~ Bi / N B^'^^ (see Supplementary Information). 
Since B'^'^'^ was several 100 mT with Bi < I mT, these 
arguments imply that the nuclear fleld distribution was 
narrowed by more than a factor of 10. Future experi- 
ments will aim at an explicit demonstration of narrowing 
via Ramsey-style experiments. 

Narrowing of the nuclear fleld distribution would 
greatly enhance our level of control of the electron spin 
dynamics. Furthermore, the observed locking effect al- 
lows us to accurately set the spin resonance frequency of 
an electron in a quantum dot to a value determined only 
by the externally controlled excitation frequency. Finally, 
our measurements suggest that we can selectively control 
the ESR frequency in one of the dots, which could be ex- 
ploited for independent addressing of electron spins in 
quantum dots that are less than 100 nm apart. 
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TUNING THE DOUBLE DOT 



The conditions for observing a pronounced electron-nuclear feedback are as follows. Qualitatively, the interdot 
tunnel coupling and the tunnel coupling to the outgoing lead are increased compared to the regime of Ref. 
Furthermore, the potentials of the double dot are tuned such that the interdot transition occurs without energy 
loss: at low power, the configuration of the dot potentials is such that electrons can tunnel elastically from the left 
to the right dot when spin blockade is lifted. Thereby, the interdot transition is made from the (1, 1) singlet to the 
(0,2) singlet, where {m,n) represent the effective electron numbers on the two dots. This working point cannot be 
used at strong driving, since the electric field component of the excitation causes photon assisted tunneling to the 
(0, 2) triplet, thereby lifting spin blockade irrespective of the spin states of the two electrons. Instead, the double dot 
must be tuned such that the (0,2) singlet electrochemical potential is higher than that of the (1, 1) singlet. This is 
nominally in the Coulomb blockade regime, but photon-assisted tunneling now provides the missing energy in order 
to make the transition from the (1, 1) to the (0, 2) singlet. 



In this section we derive an estimate for the typical magnitude of nuclear field fiuctuations around a stable point 
close to resonance. For the sake of argument we show here the derivation for a single quantum dot, although a similar 
argument holds for our double dot setup. In the double dot case, a two dimensional Fokker-Planck equation must 
be considered, where stable points correspond to zeros of {dtXi, 81X2} in the plane {xi,X2)- The results however, are 
qualitatively the same as in the single dot case. 

We consider all possible configurations of the nuclear spin system in the dot as discrete points, labeled n, defining 
n = ^{N^ — where ^+(-) denotes the number of nuclei with spin up(down) [2]. This results in iV = iV_|_ -I- iV_ 
possible values for n, ranging from —N/2 to N/2. To investigate the stochastic properties we derive a Fokker-Planck 
equation for the probability distribution function V(n), starting from a simple master equation 



In this equation V(n) gives the chance of finding the system in state n, and T±(n) is the rate at which the spin bath 
flips from the configuration n to n ± 1. We go over to the continuous limit, justified by the large number of nuclei 
iV ^ 10^ 's'l, and expand all functions around n up to second order. We find 



a Fokker-Planck equation where all rates T± are still functions of n. Due to the large number of nuclei, the spin 
fiip rates r± do not change on their full scale when increasing n by only ±1 (the features of r± occur on the scale 
of the width of the resonance ~ 1 mT, whereas changing n by ±1 corresponds to A/N ~ 5 /J.T). This implies that 
|9„r±| ^ r±, which allows us to neglect one of the cross terms resulting from the last term in (jl]). 

In the resulting continuity equation, the right-hand side corresponds to the derivative of a probability flux. In 
equilibrium this probability flux must vanish, which enables us to write down a general equilibrium solution of 
In terms of the bath polarization x = 2n/N this solution reads 



SUPPRESSION OF FLUCTUATIONS 



dVjn) 
dt 



V[n)[T+[n) 4- T^{n)] + V{n - l)T+{n - 1) + V{n + l)r_(n + 1). 



(3) 




(4) 




(5) 



Maxima and minima of this distribution are found at the zeros of the derivative of the exponent. Suppose the point 
is one of these solutions corresponding to a maximum of P(a;) (i.e. the second derivative in the point xq is negative). 
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We then expand the exponent of V{x) up to second order around the maximum, giving a Gaussian approximation for 

V{x), 



V{x) 



exp 




■dx' 



N d T, 



r+ + r_ ' 2 dxT+ + r. 

where a gives the width of the distribution. So we find that 
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d r+-r_ 



dx 



(7) 



We now only stiU want to translate this expression in terms of the 'pumping curve'. We use the relation dx/dt 
(2/iV)(r+ — r_) and define 7(0;) = {2/N){Y'+ + r_). In the limit of small polarizations, i.e. ^ 1, we can write 



^ = L{x)--i{x)x. 



(8) 



In this notation the effect of Tp (main text) is separated into two parts: (i) a polarization-dependent net spin pumping 
contribution, L{x), and (ii) a polarization-dependent contribution to the relaxation, which together with the intrinsic 
relaxation rate 1/t„ is written as j{x). One can rewrite equation ^ in terms of dx/dt and 7(0;) using the relations 
given above. This gives us finally the expression 
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(9) 



To get an idea of the magnitude of this variance, we approximate the derivative of the pumping curve at the stable 
point as roughly the height of L{x) over the width (see Fig. 4 in the main text), i.e. —dx{dx/dt)\xo ~ U'^^^/x, 
where x is the width of L{x). From equation ^ we see that we can write for the absolute maximum of achievable 
polarization x"^'^^ = L™'^'^/^{x'^'^^). Combining these two expressions and using that 7(0;™'^'^) ~ 7(2:0), we find the 
order of magnitude of the variance to be 
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N X" 



(10) 



In terms of the effective nuclear field _Bjv 1 this variance reads 



Bi 



(11) 



where ft = A/g^.BvN are the diffusive fluctuations around the unpolarized state, and Bi is the scale of the width of 
the pumping term L, in our case given by the strength of the microwave driving field. 



STATISTICS OF SWITCHING 

Here we explain how we calculated the purple curve in Fig. 4e in the main text. We suggest that the second current 
jump (red diamonds in the Figure) corresponds to the resonance being lost in one of the two dots. This occurs when 
the effective barrier between the polarized and unpolarized states becomes small enough for a typical nuclear field 
fluctuation to overcome. If we assume a simple linear decrease of this effective barrier for increasing Bf^ and include 
the effect of the finite sweep rate i?o, we find the polarization-dependent switching rate 
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where B^'^^ is the maximal sweep rate to observe any locking at all. From this expression we can derive the standard 
deviation in B^ where the second jump is observed, CTsw, and the average switching field (B^). Explicitly, we find 
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FIG. 6: Current levels of the zero-field peak and the two plateaus. An ofi'set is substracted from all current levels given by the 
average current between the tail of zero-field peak and the ESR resonance. The height of the zero-field peak is determined by 
averaging 3 points around the position of its maximum, which is determined by first averaging 10 consecutive measurements 
and determining the maximum current in the averaged trace. The current levels of the first plateau are obtained by averaging 
individual traces between the magnetic field values where the first step occurs (indicated by green circles in Fig. 4a in the 
main text) and the field where the second step occurs or the field value where the current drops to zero, if that occurs before 
the second step (red diamonds and black crosses in Fig. 4a in the main text). We require these magnetic field intervals to be 
longer than 10 measurement points (corresponding to 20mT) in order not to be omitted. The height of the second plateau is 
determined in a similar way but now by averaging between the magetic field values where the second step occurs and the field 
where the current drops to zero. The resulting heights of the zero-field peak, first and second current plateaus are represented 
here by respectively blue crosses, green circles and red diamonds for different excitation powers. 



We analyzed the set of red diamonds in Fig. 4e in the main text. From (|13p we expect CTsw to be constant in first 
approximation, which is indeed observed for lower sweep rates (100-400 mT/min). The decrease of ctsw for sweep 
rates above 400 mT/min could be a consequence of the average switching field lying too close to the resonance 
condition. Therefore we averaged the standard deviation over the first four values to find CTsw — 39 mT. Using this 
value for the standard deviation, we fitted equation (fT5|) to the data in Fig. 4e. This resulted in the fitting parameters 
^max _ 289.6 mT, B™^^ = 920.7 mT/min and 7 = 6.946 • 10~^ s~^, giving a sample correlation coefficient of 
R = 0.948; The resulting fitting curve is plotted in purple in Fig. 4e. Another way to estimate Bf}'^^ and B^'^^ is to 
extrapolate the set of red diamonds in Fig. 4e to the two axes. In this way one finds the estimates B'^^^ w 300 mT 
and B^^^ w 900 mT/min, both in reasonable agreement with the results of the ffi. 



ANALYSIS OF ESR CURRENT LEVELS 



Next to the position of the current jumps, we also analyzed the height of the current plateaus between the jumps 
as function of driving amplitude. For different microwave powers we repeatedly swept the external magnetic field 
from low to high with a sweep speed of Bq = 400 mT/min, keeping the driving frequency fixed at / = 400 MHz. For 
each trace we averaged the current of the first plateau and the current of the second plateau, and we determined the 
height of the zero-field peak. The result is plotted in Figure IH] as a scatter plot for the different microwave powers. 
We clearly observe that in all traces the highest current was measured in the zero-field peak, and that the second 
plateau exhibited higher current than the first. As to the dependence of the current levels on driving power, we see 
that (i) the height of the zero-field peak tends to decrease with increasing excitation power and (ii) the height of the 
ESR current plateaus seems nearly constant. As we attribute the observed double step feature to dragging of the 
nuclear field, first in two and then only in one dot, we here give some general considerations concerning the current 
levels during resonant electron transport in double quantum dot ESR experiments. 

Let us first consider the limit of strong microwave driving with a saturated ESR, i.e. guBBi/h much larger than all re- 
laxation and decay rates. If both dots are exactly on resonance, the driving causes the electrons to evolve entirely within 
the triplet subspace Q, i.e. in the cycle \r+) ^ {\r+) + v^fTo) + \r^)} fT^) ^ " V^ITq) + [T-)} ^ 

fr_|_). As all three (1, 1) triplet states are Pauli spin blockaded, current can only flow to the extent there is relaxation 
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from the triplets to the singlet. If only one of the two dots is on resonance, the system will evolve due to driving in 
the cycle \r±) {[To) ± p)} ^ f^i)) where in the course of every cycle the state {fro ) ± p)} can decay via 

the (0, 2) singlet to the outgoing lead, giving rise to a current. Therefore, we expect in this limit of strong driving to 
observe the highest current when only one dot is on resonance. Since the resonance is saturated in the strong driving 
regime, we expect to first approximation no dependence of the current on microwave power. 

In the limit of very weak driving, with gjjLBB\/h much smaller than the relevant rates, one would expect quite the 
opposite. In this case the system spends most time in a Pauli spin blockaded state. The blockade can be lifted by spin 
relaxation in one of the two dots or by a spin flip in either of the dots caused by the driving field Bi . In this limit we 
therefore expect increasing current with increasing driving power, and furthermore that current will be highest when 
both dots are on resonance, simply because more spin flips take place. 

During a field or frequency sweep, it is in principle possible that a nuclear field builds up in only one dot when 
the nominal ESR condition is first reached, subsequently locking the dot to the ESR condition. However, it is very 
unlikely that a nuclear field would build up in the other dot at a later time, when the ESR frequency in that dot 
is very far away from the driving frequency. A much more likely scenario is that a nuclear field builds up in both 
dots when the ESR condition is first reached (first plateau, low current), and that at the second current jump, the 
polarization in one dot relaxes to zero and only the other dot polarizes further (second plateau, high current). This 
would suggest that our experiments were performed in the regime of strong driving. 

However, there is an issue which does not fit in this simple picture. The decrease of the zero-field peak height 
for increasing power suggests that the electric field component of the excitation smears out the current peak in gate 
voltage space due to photon-assisted tunneling (at high frequencies, discrete sidebands are visible, at low frequencies, 
the sidebands overlap). This could account for the decrease of the zero-field current, but should presumably aff'ect the 
ESR current levels in the same way since the ESR transition is saturated at strong driving. However, experimentally, 
the current levels at the two ESR plateaus are roughly independent of power, rather than decreasing with power. 
This point remains at present unresolved. 

In order to develop a coherent picture of electron transport at zero-field and at spin resonance, a more systematic 
and detailed study of the dependence of the current levels on driving power and on the tuning of the double dot 
(tunnel coupling, detuning) is needed. This is quite involved, since the behavior of even the zero-field peak varies 
widely with tuning parameters. 
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